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Key Takeaways 


e New diffuse equations derived from same 
assumptions as GGX-- Smith specular 


e New cheap and good shadowing/masking 
function G for GGX-4 Smith specular 


e New trig identities for shader optimization 


| SA 


Some fun discoveries on the way 


e Why's that “4” in PBR specular? 

e What is “s” in Oren Nayar diffuse? 

e Smith shadowing/masking assumptions 
e A physical interpretation of Lambert 

e Help interpreting Disney’s BRDF slices 


à ^ 
Quick aside: Original motivation 


e Titanfall 2 used Oren-Nayar diffuse 
e Question: How to get Oren-Nayar's 
roughness s from GGX's roughness a? 


e Discovery: Oren-Nayar came from very 
different assumptions! 
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Quick aside: Original motivation 


Shadowing/Masking V-cavities Smith 
Normal Distribution Spherical Gaussian GGX 
Roughness parameter s € [0,00] a € [0,1] 
Perfectly flat s0 a=0 
Standard deviation of S 0 a=0 


2 
Slopes of normals a^o a#0 
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Quick aside: Original motivation 


e Oren-Nayar and Smith+GGX don't match! 
e Can't even match standard deviations 
e Hmm... GGX standard deviation is aoo 
e Maybe “best” to mipmap/filter a?? 


e Sum of two GGX distributions is not GGX, so can't 
mipmap/filter "properly" 
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Road map for today's talk 


e General microfacet-based BRDFs 


e Simulating diffuse for GGX-4 Smith 
microfacet model 


e Comparing to other diffuse BRDFs 


| FEB 27-MAR 017 | EXPO: MAR 1-3, 2017. #GDC17 


Road map for today's talk 


e General microfacet-based BRDFs 


e Simulating diffuse for GGX-c Smith 
microfacet model 


e Comparing to other diffuse BRDFs 
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Microfacet BRDF sub-topic map 


e General form 
e How we get PBR specular from that 
e Extend to diffuse BRDF 
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Microfacet BRDF sub-topic map 


e General form 
e How we get PBR specular from that 
e Extend to diffuse BRDF 


à SA 


Microfacet models 


e Complex macrosurface BRDF averages 
many that use a simple BRDF 


e Basically just subpixel shader antialiasing 
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Real world examples 


Images: http://funjungle.net/the-world-is-different-under-the-microscope/ 
http://www.geek.com/news/chocolate-under-an-electron-microscope-looks-like-an-alien-planet-1648301/ 
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"General" Microfacet-based BRDF 


e Not fully general; assumes heightfield 
e No weaves, arches, or caves 
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General Microfacet-based BRDF 


(m-L) (m-V) 


o fo Pm(L,V,m)D(m)G2(L, V, m) == 


IN-L| |N-V| 


General Microfacet-based BRDF 


(m-L) (m-V) 
IN-L| |N-V| 


e f, Pm(L,V,m)D(m)G2(L, V, m) == 


e Integral over all microfacet normals 
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General Microfacet-based BRDF 


(mL) (m-V) 

e Jo D(m)G>(L,V, a rar N V] 

e How an individual facet , 
responds to light 


e I.e., microfacet BRDF 
centered on m 

e Usually ideal mirror or 
ideal diffuse 
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General Microfacet-based BRDF 


(m-L) (m-V) V) 
Jo Pm(L, V, m) G2(L, V, m)-—— IN-L| |N- VI 
e Probability density of normal m d 


e Which facet normals exist, but 
not their arrangement (shape) 
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General Microfacet-based BRDF 


e Jo Pm(L,V,m)D(rm) fear Im 

e Occlusion due to actual microfacet & 
arrangement (actual shape) 

e A.k.a. shadowing/masking function 2 cand 

e Probability microfacet m sees both Ss. 

light L and viewer V Ee 


Microfacet BRDF: G, vs G4 


e G;(L,V,m) is 9o visible in 2 directions 
e G (V,m) is 9o visible in just 1 direction 
e In practice, G» is derived from G, 


G4(V,m) 
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Microfacet BRDF: D(m) properties 


(m-L) (m-V) V) 


9 Jo Pm(L, V,m) G2(L, V, m IN-L| |N-V| 


e Probability density of normal m 
e How quickly cumulative probability changes near m 
e Will change more quickly in more probable regions 
e In range [0,00], not [0,1]! 
e D(m) = © for any m whose probability + 0! 
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Microfacet BRDF: D(m) properties 


e Total surface area of all microfacets 
e Always > 1 if any roughness at all! 
e Jo Damn) cos O, dm = 1 


e To normalize total area, project microfacets 
onto macrosurface using cos0,, 2 m: N 
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General Microfacet-based BRDF 


(m-L) (m-V) 
è Jo Pm, V, m) IN-L| |N-V| 


e Probability density of having microfacet 
normal mthat is both lit and seen 


e I.e., probability density of using BRDF pm(L, V,m). 
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General Microfacet-based BRDF 


(m-V) 
© Jo om (5, V, m) D(m)G; CL, V, m) TT 
° - How big facet m appears to the light 
o El - How big facet m appears to the viewer 


e I.e., normalize contribution and to viewer 
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General Microfacet-based BRDF 


” Jo Pm(L, V,m) dm 


e Probability density of light from L reaching V 
in a single bounce off microfacet normal m 


e Requirement: f, dm € 1 
e Only = 1 for flat D(m) - too dark if any roughness! 


ws 


General Microfacet-based BRDF 


e Related requirement: 
| vena, m)&m - V)dm = |N -V| 
Q 


e In any direction V, total visible microfacet 
area equals macrosurface area 
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Microfacet BRDF sub-topic map 


e General form 


e How we get PBR specular from that 
F(L, H)D(H)G4(L, V, H) 
4IN - LIIN- V| 
e Also: Where does that 4 come from, and why isn't it z? 


e Extend to diffuse BRDF 


»p 
PBR Specular Microfacet BRDF 


e Microfacet BRDF is a perfect mirror 


e I.e., light reflects if and only if m = H 
e Mathematically, BRDF is a scaled dirac delta ó,,(H,m) 


X leas Leah 


| o 4h 
/ sê ER) 
~ \Y UBM 


PBR — Microfacet BRDF 


e Pure mirror BRDF: kó,,(H,m) 
e ôm(H,m) is the dirac delta using measure m 
e k is some normalization factor we must find 
e Normalized BRDF: J, o(L, V, N) cos6y dV = 1 


e For any light and normal, all energy reflects 
to exactly one viewer 


Specular BRDF normalization 


e General case: J, p(L,V,N) cos0y dV — 1 


e Our case: J, kó,, (H,m) cos Oy = 


e Must integrate over dm to evaluate 6,,, so find 
dV í : . 
am tO change integration domain 


Specular BRDF normalization 


e — is how fast V changes relative to m 


e This will introduce PBR specular's 4! 
e Next few slides show how 


Specular BRDF normalization 


e We're going to find dm from dV to get T 
e First, 6,, picks m = H, SO dm = dH 


e All vectors sketched 
on unit hemisphere dV 


Specular BRDF normalization 


e Move solid angle dV... 


dV 


Specular BRDF normalization 


e Move solid angle dV to L +V... 


Specular BRDF normalization 


e Move solid angle dV to L+V, scale by... 
e L+V sphere: H-V 


Specular BRDF normalization 


e Move solid angle dV to L +V, scale by... 
e L+V sphere: H-V . 
DI a Do sil 


e Unit sphere: en ee 


Specular BRDF normalization 


e Move solid angle dV to L+V, scale by... 
e L+V sphere: H-V sav 


: , 4m 1 
e Unit sphere: a eae 


H-V 


IL+V|2 av 


e dm = 


Specular BRDF normalization 


e Move solid angle dV to L+V, scale by... 
e L+V sphere: H-V 


: , 4m 1 
e Unit sphere: a eae 


H-V 


IL+V|2 av 


e dm = 


Specular BRDF normalization 


eIL+VI=H-(L+HV)=H-L+H-V=2H-V 
e This 2 squared is specular BRDF's 4! 
H-V H-V 


e = —— = ————— L + V 
dm Lave dV EV? dV ; 
am 


Specular BRDF normalization 


e. Jokôm(H,m) cos 0, — dm = 1 


° I5 kó4,(H,m)(m -V)(4H -V)dm = 1 
1 


————— since m = H and H-V-—H-L 
A(H-L)(H-V) 


ek= 


Om(H,m) 
4(H-L)(H-V) 


e So, pure mirror BRDF: 


à SA 
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Specular Microfacet BRDF 


17 | EXPO: MAR 1-3,2 


e Only the Fresnel reflection fraction F(L,m) 
of incoming light does specular reflection 


e So, final specular microfacet BRDF: 
ôy (m) 
IH - L||H -V| 


pm, V, m) = F(L,m) 
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Specular Microfacet BRDF 


(m-L) (m-V) 
© jo Pm Us V, m)D(m)Ga, V, m) S NV 
F(L,m)6m(H.m) (m-L) (m-V) 
Jo 4|H-L||H-V | dm 


e ôm(H,m) eliminates integral and sets m = H 


e Specular BRDF: 


F(L,H) 
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Microfacet BRDF sub-topic map 


e General form 
e How we get PBR specular from that 
e Extend to diffuse BRDF 
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Diffuse Microfacet BRDF 


(m:-L) (m-V) 


© fo Pm(L,V,m)D(m)G2(L, V, m) —— 


e Lambertian diffuse: o,,(L, V, m) = - 


IN-L| = VI 


e No dirac delta to eliminate integral & 
e No closed form solution for GGX-« Smith 66 
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Diffuse Microfacet BRDF 


e Solved integral numerically, hoping to 
find good approximation 
e Same approach as the Oren-Nayar paper 
e Up to half the light was missing! 


e Can't ignore multiple bounces... 
e (Full Oren-Nayar includes a second bounce too) 


Importance of 


à LE E > e. - proper albedo 
| | Top: 
"a ` ` . ) , Correct 


| Bottom: 
135º DO | i } | Albedo x 0.5, 
Light x 2 


^ am og 
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a= 0? 0.252 0.52 0.752 ile 


Importance of 
proper albedo 
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| Top: 
"a ` ` > ) ) Correct 


| Bottom: 
=s E | i | | Albedo x 2, 
Light x 0.5 


^ am og 


a = 0? 0.252 0.52 0.752 ile 
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Road map for today's talk 


e Simulating diffuse for GGX-« Smith 
microfacet model 
e Comparing to other diffuse BRDFs 
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Road map for today's talk 


e General microfacet-based BRDFs 
e Simulating diffuse for GGX+Smith 
microfacet model 
e Shadowing/masking functions 
e Path tracing 
e Comparing to other diffuse BRDFs 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 
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Height-Correlated G 


e G, is geometric visibility to n directions 
e If uncorrelated: 
G;,(L,V,m) = G4(L,m)G4(V, m) 
e Not realistic! Higher points more likely visible 
to both L and V (and lower points less likely) 


e Still, surprisingly good in practice 


Height-Correlated G 


e Uncorrelated G takes light 
hitting a normal in the 
heightfield... 


S TUM 


zd 


93% 


8790 
57% 
0% 
0% 


Height-Correlated G 


e ...and redistributes it evenly 
across each microfacet with 


that normal 


NT 


zt 


93% 


8790 
57% 
0% 
0% 


,2017 | E :M 


Height-Correlated G 


e This tends to move light 
lower, reducing its visibility -2 9396 
and darkening specular. -1 | 87% 


Height-Correlated G 


,2017 | E :M 


à SA 


e Uncorrelated G's error is related to 
occlusion 


e Error bigger for rougher surfaces 


e Error bigger when L and V more glancing 
e No error if L= N and/or V =N 


VT 
: i : | 
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Uncorrelated vs. Correlated G 


Height-correlation 


’ | 135° (bottom) boosts 
glancing reflection 
| | | — ass On rough surfaces 
| | Tr B Black albedo; 
light intensity = 7 
(DUE 0.25? 0.5? 0.75? 12 
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Uncorrelated vs. Correlated G 


G ph | 
| ^^ q | 
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Exact vs. Approx Correlated G 


Approximation 

ww Ww ss (bottom) is quite 
good, but still a 

| | | 157.5? little too dark for 
medium angles 

| ) | ,, and roughness 


eS 0.25? 0.52 mms E 
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Exact vs. Approx Correlated G 
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Height Correlated G Uncorrelated G Difference 


dd 

G 

e There is angular correlation too 
e L=YV should have: G;(V,V,m) = G4(V,m) 
e Uncorrelated form: G,(V,V,m) = G4(V,m)* 


e Height correlated G; somewhere in 
between when L=V 


Correlated 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 


Smith Shadowinia/Maskina 


e Assumes all normals equally occluded 
e I.e., G and G, don't depend on m 
e Most balanced assumption possible 


ME DEVE ER 


Smith Shadowing/Masking 


e Can derive from normalization constraint: 
6,Q) | Dmm - V)dm = |N -V| 
Q 


à SA 


e Can also derive from ray-tracing a 
probabilistic heightfield 


GOC GAME DEVELOPERS CONFERENCE | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017 #GDC17 ~ SA 


Smith Shadowing/Masking 


e Super basic ray trace derivation: 
e Project m onto 2D plane with PDF P5;(p, q) 
e D(m) isotropic, so use 1D slice with PDF P,(q) 


e Use P,(q) to get PDF of ray-surface collisions 
while ray with slope u walks the heightfield 


e Use PDF of collisions to get G, 


| FE 


Smith Shadowin 


g/Masking 


N 


Polar normal m; PDF D(m) 


| FE 


Smith Shadowin 


pr Ye 


g/Masking 


2D slope p,q; PDF P>2(p, q) 


N 


Polar normal m; PDF D(m) 


| FE 


Smith Shadowin 


g/Masking 


2D slope p,q; PDF P>2(p, q) 


q 
(p, q, Ym 


All p/for q Y 


Polar normal m; PDF D(m) 


Smith Shadowing/Masking 


2D slope p,q; PDF P55(p,q) 


q 1D slope q; PDF P,(q) 
(p,q, 39e 


All p/for q v 


Polar normal m; PDF D(m) 


Smith: Arbitrary D(m) 
e P22(p,q) = cost Om D(m) 

© P2(q) = f. P2; p, q)dp 

e Au) = AC — H)P2(q) dq 


e GV) = 


1 
1+A(L)+A(V) 


Go(L,V) = 


EC 
1+A(V) 


Smith: Correlated vs Uncorrelated 


1 


e G(V) = 1+A(V) 
l IEEE XN 
e Correlated: G2(L, V) = —— OA 
l COME QEMIME 
e Uncorrelated a G> (L, V) m 1+A(L)+A(V)+A(L)A(V) 


e Too small, unless A= 0 (i.e. G, = 1) for L or V 


Smith for GGX: e 


| Do ato 

e For GGX: D(m) = T (cos* On(la2+tan? 0,4)?) 
2 A 

© Pop) =——— = : 


n(a?-tan? 0m)?  n(al+p2+q?) 
2 


a a 
* PQ) = farra P = Marea 


| FE 


Smith for GGX: A(V) 


1 [oo 1 2 +u? 
+ AQ) = f (a Pada EE (3 - 1) 
e u = cot0y 
e cos0y =N-V 


o ACV) = 1(38 8-020 PP q 
2 


N-V 


Smith for GGX: "Tr "7 V) 


2N.V 
A a2 +(1—a2)(N-V)2+N-V 
. G(L,V) = 2(N-L)(N-V) 


N-V4/a?--(1—a22)(N-L)? -N-.LJ a? - (1-a?)(N-V)* 
e Would like cheaper approximation! 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 


R 3, 2017 | EX 


mate GGX G,(V) 


e Denominator of G,: 
e Ja? -(1—a?)(N-V)? - NV 


° lerp((N -V)2,1,a2)+N-V 


e Approximation: lerp(N-V, La) +N-V 


Smith: Approxi 


017 | EX 


Smith: Approximate GGX G4(V) 


«GO = 2N-V | 2NV 
1 lerp(N-V,1,@)+N-V N-VQ-a)+a 


e Turns out, identical to Unreal's Smith: 


N-V a 
e GV) = EE 2 
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Smith: Approximate GGX G,(L,V) 


e Solve this G, for A(V), plug in for G;(L,V): 
2|N-L||N-V| 

lerp(2|N-L||N-V |,|N-L| - | N-V],a)) 

e G, S numerator cancels in full specular BRDF: 


e Go(L, V) = 


4|N-LIIN-V| 2 lerp(2|N-L||N-V],|N-L|+|N-V],@) 
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Smith Approximation Cost 


e Compare cost of denominator: 
F(L,H)D(H) 
(IN-L|(2-00)+0)(IN-V|(2-0)+a) 
F(L,H)D(H) 


e G4(L)G4(V): e 4 cycles 


e G;(L,V): ~6 cycles 


e Costs exclude calculating N-L and N-V 
e Height-correlated form has negligible extra cost 
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Smith Approximation Quality 


e Helps rough dielectrics at glancing angles 
Exact G2(L,V) Approx G;(LjV) Approx G, (L)G,(V) 


Top left of image: 


Fo = 0,1 = 20 glancing reflection 


Bottom right: 
Ry elle2 normal incidence 
and viewer 


GGX Specular BRDF for a = 0.8; N- L,N-V increase down,right 
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Smith Approximation Quality 


e Helps rough dielectrics at glancing angles 


Exact G2(L,V) Approx G;(L,V) Approx G1CL)G4 (V) 


” glancing reflection 
Bottom right: 

F=? normal incidence 
and viewer 


GGX Specular BRDF for a = 0.8; N- L,N-V increase down,right 


Top left of image: 
Fg m 0,1 = 20 


—- yar 
se Be uat 
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Uncorrelated G 
Approximation 


É cd e mi 


Difference Image: 


Red - correlation 
Green = approximation 
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a 


EC 


Fe 


Correlated G 
Approximation 


Difference Image: 


eS 


Relative to 
uncorrelated 
approximation 
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——. GEES IEEE 
= sé Exact 


Difference Image: 


Relative to 


correlated 
ras approximation 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 


ME DEVE ERS 


Smith Microsurface Heightfields 


e Example 1D heightfield 
from Smith ray tracing 


derivation sketch in 
Walter 2007 


| SA 


Smith Microsurface Heightfields 


e Derivation *also* uses 1D heightfield of 
mostly independent slabs 
nearing zero width T> 


e Only forbids suddenly 
being under heightfield 
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Why Smith masking is weird 


e Ray tracing derivation has contradictory 
assumptions at different steps: 


e Height in next dr independent of this height 
e Assumes not continuous 


e Heightfield is any differentiable function 
e Assumes continuous 


coo ooo DA Qe ^ 
Why Smith masking is weird 


ER E | FE 2 | EXPO: MAR 1-3, 201 


e Math says visibility is asymmetric: 
downward rays less likely than upward 
rays to survive the same heightfield path! 


e AQ) = 7 f, (4 — u)Po(a) da 


e A(u) integrates all q > u, so u < 0 can hit more 
values of q than when u » 0 
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Why Smith masking is great 


e Only energy conserving G where all facet 
normals have the same fraction visible 


e Any other G not using m gets total visible 
area wrong for some directions 
e Too high reflects too much, creating energy 
e Too low reflects too little, absorbing energy 
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Diffuse Simulation sub-topic map 


e Shadowing/masking functions (G,, G;) 
e Uncorrelated vs height correlated G 
e Smith shadowing/masking 
e New Smith+GGX G, approximation 
e Greatness and weirdness of Smith 


e Path tracing 


à SA 


Path traced diffuse solution 


e Smith oddities prevent real heightfields 
from matching its assumptions 
e Can't be both continuous and discontuous 


e Must ray trace the mathematical model 
e See bonus slides for numerous details 


à SA 


First ray traced result 


e Simple ray tracer with Fresnel to choose 
GGX specular or Lambertian diffuse 

e Resulting BRDF was not symmetric! 
e p(L,V, N) + p(V,L,N) 

e What went wrong? Both parts are 
symmetric BRDFs! 


Cause of asymmetric BRDF 


e Essentially had a merged BRDF: 
e p = F(L,N)Pspec + (1 — FC, N))pairy 
e Fresnel interpolation asymmetric! 
e How to fix in a physically plausible way? 


GOC GAME DEVELOPERS CONFERENCE | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017 #GDC17 | 4/ 


Why Lambertian diffuse? 


e What does Lambertian diffuse simulate? 
e BRDF p = =: same for all viewers 


e Radiance = pcos 86y: more photons at normal 
1 
cos Oy 


e Balanced by for total surface area seen by V 


e Why the cosine energy falloff? Answer is surprisingly 
hard to discover, yet quite simple! 


GOC GAME DEVELOPERS CONFERENCE | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017 #GDC17 On ^ 


Lambertian diffuse explained 


e BRDFs given at the surface, 
but diffuse light just passes 
through the surface 

e Lambert assumes interior light 


has same density in all directions 
e Cosine falloff is from surface ce Ae > > 
angle relative to light direction Xe Xe e x 


" "d 
N 
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Lambertian diffuse explained 


e Same energy per area each 
direction 
e Directions angled to surface 
project over larger area 
e Area per unit light scaled by !/.,. 4 
e Light per unit area scaled by cos 0 
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Lambertian diffuse explained 


e Light enters, bounces 
around, exits 

e Exit direction is random 
after many bounce events 

e Albedo comes from 
frequency-dependent 
absorption events 


| SA 


r symmetric BRDF 
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Fixing diffuse fo 
e Entering uses Fresnel for 


reflection/transmission 


e Exiting assumes always 
transmit 


e Exiting needs Fresnel too! 


| SA 


r symmetric BRDF 
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Fixing diffuse fo 


e Reflect: F = F;-(1— F(1—-N-V)? 
e Transmit: 1— F 2(1— Fj(1—-(1—N-V)?) 
e Fresnel's laws are symmetric, so fraction 


entering surface from viewer equals 
fraction exiting surface toward viewer 
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Fixing diffuse for symmetric BRDF 
e Internally reflected light keeps getting 
chances to transmit; need to normalize! 


e 21 pr” k(1 — (1 — cos @)°) sin dð = 1 


e Factor (1 — F9) absorbed into norm factor k 
o needed to normalize a BRDF 
e 27 and sin 0 dé from integrating on a hemisphere 
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Fixing diffuse for symmetric BRDF 

e Easily solved exactly: k = — = — 

e Merged diffuse+spec microfacet BRDF: 
eF=F,+(1—-F))(1—-N-V)° 


1.05 
ki P =F Pspec + (1 F) —(- (1- NV?) 


Finally! 


e Now have everything needed for path 
tracing simulation, resulting in... 


A 
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| Approximate 
. * *. diffuse only 


0° 
" ` ) ) (0.75,0.5,0.25] 
135° 


a=07 0.252 0.52 0.752 ile 


VEO am Cs 
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GGX Diffuse Approximation 


facing = 0.5+4+0.5L-V 
rough = facing(0.9 — 0.4facing) e) 
smooth = 1.05(1— (1— N- L5) -(A- N -V)9) 


single — = lerp(smooth, rough, a) 


multi = 0.1159« 
diffuse = albedo * (single + albedo * multi) 


Aside: Useful shader identities 


e |L4+V\2=2+2L-V 
i 0.5 + 0.5L: V =Ż|L VP? 


e N-H= N-L+N-V 


IL+V| 


; L-H=V-H=-|L+V| 
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Aside: Useful shader identities 


e Can find N -H and L- H without finding H! 


Get H = normalize(L + V) 
Get H then N-H 
Get H then N- H and L- H 
N - H from identities 
L:H and V -H from identities 
* Add 3 cycles if you don't already have L- V 


#GDC17 


13 
16 
19 
7x 
g* 


NN A A AȘA 


Aside: Useful shader identities 


e lenSq LV 22 2L: V 

e rcpLen LV = rsqrt(lenSq LV) 

e N-H-(N-:L«-N-:V)xrcpLen LV 

e L-H =V -H = rcpLen;y t rcpLenjy * L-V 


: 1 1 _ 1 24+2LV\ — l 
; (Since L: H - 2|. V| 2242 €21-Y =EL) = (1+1 V) 


1 
mu 
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Road map for today's talk 


e General microfacet-based BRDFs 
e Simulating diffuse for GGX-c Smith 
microfacet model 
e Shadowing/masking functions 
e Path tracing 
e Comparing to other diffuse BRDFs 


But First... 


e It's good to quickly understand Disney's 
BRDF slices 


à SA 
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Disney's BRDF slices 


e BRDF is a 4D function of 2 polar vectors 
e Before, light+viewer vectors: 0,, 4,,0,, dy 


e After, half angle+difference: 6,, dp, 0a, Pa 
e Isotropic BRDFs never depend on q, 


e Dependence on 4$, is often negligible 


à SA 


Disney's BRDF slices intuition 


e Each row is a light+viewer pair (0,) 


e Opposite at top, perpendicular in middle, 
coincident at bottom 


e Left-to-right shows falloff going away 
from center of specular highlight (6, ) 
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False color example on lit sphere 


BRDF Slice Corresponding Lit Spheres 


Lighter bands highlight Lighter bands highlight bp; = 90º; 
rows used by spheres ba increases counterclockwise 


Disney's BRDF slices 


Silhouette Opposite 


Perpendicular 


‘ Retroreflection 


Equal 


Behavior of 8n, 804, pa on spheres 


ONE 


0,;N-H 05 L-H;L-V 


Disney's BRDF slices 


e Various identities: 
e cosÜü, — N-H 
e cos; =L-H=V-H cos204 =L-V 


da = N-V-N-L 
° COS Pa = J(2-2L-V)(A-(N-H)2) 
e BRDFs mostly functions of N- H and L.V! 


Almost ready to compare BRDFs! 


e First, introduce the comparison format 


is 
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Title says which diffuse 
model is shown. This 
intro uses the new model 


V mm O: 
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) o) 4 ) À < This panel shows the 
same lit spheres as 
previous examples. 


y 
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OCLI. New Model 
L229) 


The matching BRDF 
— slices are here 
<— (uncorrelated G) 


^— a oa 


. PAN 
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e^o 


monem 
^Y Oy 03 ^ 


VO m O: 


uffizi 


Same full BRDF with « 
from O to 1, but lit by 
Paul Debevec's HDR 
environment probes. 


campus grace 
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s. e o e e Lambert 
€COC nossa 
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a ¥ T T F 
à r 
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y n o . . New model 
$t | AMAM 


p 
"— 
A " 
l | 
| 
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. > e o e New model (hybrid) 
Smooth uses Disney's 


] . ` ) ) faoo = 0.5, SO same as 


Disney when a = 0 


“em: 


ED... resmas or é € 
* > e ^ o New model (cheaper) 


Smooth uses Lambert 
1 ) 


l AN 
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Lambert 


ZN 
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m 


Z N 
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New model 


Lambert 


"T 
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New model 


As 
OP o 
UBM 
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New model 


PS AY 


By. I. Ww hx 


UBM 
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Appendix 


The following is a bunch of derivations for Smith shadowing/masking 
from the ray tracing formulation, and how you use that to actually do the 
path tracing. This is how I got the results included in the preceding 
presentation. 


This is quite math heavy. As such, it fits much better in an appendix than 
in the talk. It is hard to read derivations to an audience, and it is even 
harder to listen to them! It's better to be able to go at your own pace, 
and to be able to flip back and forth as needed. 


Final caveat: I didn't polish these appendix slides much (e.g., there is a 
complete lack of figures). Still, the information and derivation should be 
helpful to those who like to understand where things come from, and/or 
who want to understand the Smith shadowing/masking derivation. 
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Solution - Path Tracing 


e Shoot photons into the microsurface for a light direction 
e See which view direction those photons come out 
e This lets us model diffuse and specular interactions 


e But first, we have to be able to ray trace the microsurface 

e The microsurface is implicitly defined by the normal 
distrubution function D(m) and the shadowing/masking function 
G(L,V,N) 

e G(L,V,N) derived from D(m) is basically ray tracing 

e SO, we need to understand how Smith G(L,V,N) works 
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Starting to derive Smith 


e This basically follows Appendix A in Walter’s 2007 GGX paper 


e With many missing details filled in, and slightly reordered. Any 
ae with Walter’s appendix are my own attempt to complete the 
erivation. 


e It may be handy to pull up Walter’s appendix as you follow these slides 
e Consider ray tracing a 2D slice of the heightfield in the plane of 
the ray and macrosurface normal 
e Y axis is height (E), X axis is projected distance along ray (t) 
e We need probability of hitting height field, given that we 
haven't yet 


e Let P (č) be the probability density of height ë 


e Probability height € is above the heightfield f(é) = f^. Pi GO dx 
e This is total probability of heightfield being lower than 2 
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PDF of hitting heightfield (1/2) 


e Forray £g t ux to hit in the next Ar from height E with slope q: 
e éotut»é 
e €& +ut+pAtT>é+qAt 
e Rearranging, we have: 
© So tut2é$»é4o ut —(q—u)N 
e Clearly requires q > u 
ipe "TT ; co régtUT 
e Probability of hitting is li SE yu t-(q-uar 1 GO P2 (Q)dx dq 
e P,(x) is probability density of height x 


e  P,(q)is probability density of slope q 
e Product assumes probability P,(x) and P,(q) are independent 
e All normals equally likely at each height; heightfield is fractal, not like canyons or spikes 
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PDF of hitting heightfield (2/2) 


co Sotut 
° Jn Pala) fe eu c- (a2 iac i ax dq 


e Can pull P,(q) out of the inner integral because it doesn't use x 
e Take dim At — dt 
T 


e Assume P,(x) is constant over dt at P,(£g + uT) 
© Ji Po(q)Pr o +u T)(q — idt dq 
Probability of ray & + ur hitting in next dr: 


dt P(£9 tut) | (q — u)P5(q)dq 
u 
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Meet S, the surviving fraction 


e Need conditional probability given we start outside the heightfield: 


dt P (o4 T) f; (a-1)P4(q)dq — dt Py(&+u T) f (a-H)P>(a)da 
° poot#® p. (ax E fGo*u v) 
e Let S(é,u,T) be the surviving fraction and consider how it changes: 


Pa (otu T) fj^ Px(q)(a-10dq 
TENS 
f Gotu T) 


e as = (de 


e Le., the fraction of surviving rays hitting in the next dr are subtracted from S 


as (se Sy orn) r 


flotu T) 
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Start solving diff. eq. for S 


as (PiGotu 0 f, Po(a)(q-n) aq 
e We have - - (Bem ems 


d - E 
e In general, "m g( = —e-9(0g'(7), soS— e IM 
rn PaGotH# 2) o Pola)la-u)da 
* 4095 fGs*u T) 
«Recall definition f(£; + u t) = f^" ^* P, (x)dx 
C = P + ut) Go tur = HPS tnr) 


«Define A(u) = J, (a — u)P5(q)dq 
eAssumes u 0 


uP: ou Da Ju Pola)la- nda 
fotu T) 


= A(u) f ote 1) 


* g()- fotu o) 
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Final solution for S(Z,, u, T) 


e g(t) = fo g'COdt = f; AQ) Testu dt = AQ) In f + u Ol 


e Uses the fact that “In g(t) = "TO KY 
e Assumes f(2) and P,(x) are continuous 


e gt) = AQu)nf(do + ux). —1nf(49)) = A(u) In p 


-AQ)In or (otro 0" 
f (Eo) 


e Slo mT) =e 9 =e 


e We've derived the heart of Smith shadowing/masking: 
S (£o, LL, T) = (Crs 


f (£9) 
f Go + ux) 
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Solving for G,(u) given S(£e, u, T) 


e We can use this to find the probability of a ray escaping 


ACu) ACu) 
© Son) = lim SGo T) = (E) = (ER) pago 


e Assumes y > 0 in f (o + ux) > f (œ) 


e Can find probability of seeing the heightfield in slope u = a. 


at 
© GQ) = f, PE) SE wg 
e Integral over all heights of the probability of having height € and escaping in 
direction u 


© GW = [^ f'G) FOAM aE = sf | 


e f (o?) = 1 and f(—oo) = 0 regardless of choice of P,(&), which defines f (ë) 
1 
1+A(u) 


« Giu) = 
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Solving for G> (ur, uy) 


e Can also find visibility in two directions from one height 
e G;(Hp Ay) = An P,(6) SG, uj SCE, uy )d$ 
© Galur uy) = So PAE) FEM f )^ dE 


e G(uL, uy) = jJ Py (E) f (p) 0t ^o gz 


e This gives the height correlated Smith shadowing/masking 
function: 
1 


e G5 (Ur, Hy) = 14 AQu 4 AQ) 


e Note that u = cot0, where 6 is the angle from the macrosurface 
normal 
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Starting to derive A(z) 


We still need to finish derivation of A(u) = "MC — u)P;(q)dq 


This uses the probability of a tangent slope P,(q), where q = E 


at 
We have the surface area of a microsurface normal, D(m). 


Project the microsurface normal area onto the macrosurface to 
get the probability density of a normal per unit area of the 
macrosurface 

e D(m)cos0,, where cos 8m =m-N is the angle from vertical 
Project from spherical coordinates (0m, bm) to plane (p,q, 1) 


(cos Pm sin Om sin Pm sin Omcosêm) _ (cos o tan 8m , sin bm tan 8m, 1) 
| m m» m m 


cos Om 
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Change variables (8m, Øm) to (p,q) 


e p = cos Ọm tan Om, q = sin n tan Om 
e Implies p? + q? = tan? 6,, 
e Probability density of normal m is D(m) cos 6, dm 
e dm = sin m dO4 dog, 
e We need it as probability density of slopes P,,(p, q)dpdq 
e This is the same as D(m) cos 0, dm, just with change of 
variables. 
e Need the Jacobian, based on partial derivatives 
Op _ a dp. _ cos dm 
di rm sin dy, tan O, Bo. GO, 
dq dq | sindm 


— COS Pm tan Om» 005, pF 0,, 


Om 


GOC GAME DEVELOPERS CONFERENCE | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017. HGDCI7 


Final Jacobian for (8m, Øm) to (p,q) 


Op Op 
Odm 290m = | ' sin $4 COS bm 
e Idet aq aq || = |— sin Pm tan Om zos? 0%, COS Ọm tan Om EIN 
dpm Om 


e This is the Jacobian for the change in area of the measure for a change 
of variables from (0,,, m) to (p,q) 


: "EN tan Om EN 2 tan Om | _ tanOg _ sinôm 
e Jacobian | sin” Om cos? Om cos" Pm cos? Om] cos? Om cos? Om 
e This means 
sin 0 
dir ddom = dpdq 


e sin6 ddom = cos? 0 dpdq 


GOC GAME DEVELOPERS CONFERENCE | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017. HGDCI7 


Completing (8,,, bm) to (p,q) 


e Change of variables for m from (6,,,¢,,) to (p,q) has 
e dm = sin 8n dO, d$, = cos? Om dpdq 
e We want P,,(p,q)dpdq = D(m) cos0,, dm 
e PD,,(p,q)dpdq = D(m) cos 8m cos? 0,, dpdq = D(m) cos* 0m dpdq 
e Recall that p? + q? = tan? 04... SO 0, is a function of (p,q) 
e If D(m) doesn't depend on m, D(m) cost 6,, dpdq is a function on 
(p, q)! 


e For GGX, D(m) = : 


T (cos^ Om(a2+tan? 0,,)?) 


2 


a? a? 


m(a2+tan? 604)?  m(a2+p2+q2)? 


e D(m) cost 6, = 
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Slope PDF from 2D to 1D 


We're getting close! We have probability density of 2D slope 
(p, q): 
e P55(p,q)dpdq 


We need the 1D probability of slope: 

e P;(q)dq 
Since we've assumed D(m) doesn't depend on m, we can 
arbitrarily rotate (p,q) such that q aligns with the ray direction 
and p is perpendicular to it 


This lets us integrate P,,(p, q)dpdq over all p to get P,(q)dq: 
© Pj(q)dq = f, P22 (p, q)dpdq 
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P,(q): normals or tangents? 


P,(q) was derived as the probability density that a microfacet 
normal goes q units along the x-axis (t) for every 1 unit along 
the y-axis (é). 

e Tangents are always perpendicular to normals. 

e In 2D, vector (x,y) is perpendicular to (—y,x) and (y, —x). 
So, this is equivalent to the microfacet tangent slope going —q 
units along the y-axis (€) for every 1 unit along the x-axis (7). 

e This means P,(q) is the probability density of tangent slope -q. 

e But D(m) doesn't depend on @,,, so P;(q) = P> (~q). 
This means P;(q) is the probability of a microfacet tangent 
slope q. 

e This is how we used it earlier 
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Use GGX's P,,(p,q) to get its A(z) 


2 
e For GGX, we saw that P55(p,q) = CEDO 

ar a 
e P, (q) = = i P22(p, q)dp = i n(a2+p2+q2)2 dp = n(a2+q2)15 


e All this is to find A(u) = "M (q — 2) P2(q)dq 


_ : o a?(q-L) 8^ *qu _ 1, ya?sg? 
" DOS le ugue - 2p a? q? gy PA És 
e We have u = z * for a view vector, so p = cot 8y = y * for 8, from 
vertical. 


1 a? sin? 0y -cos? Oy a? * (1—a7?) cos? 0y 
Au) = — À + dE sil Dy sto Oy L ye i-a) cost dy 


2 cos Oy 2 cos Oy 


1 
2 
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Smith masking is weird (1/2) 


e Ray-tracing derivation of the Smith masking function assumed 
any height/slope can be immediately adjacent to any other 
height/slope. 

e Le., the heightfield is continuous nowhere, yet differentiable 
everywhere. 
eTo be differentiable, you have to be continuous, so this is contradictory. 


e This also means id can't integrate slopes to get heights, yet slopes are 
the derivatives of the heights. 


eThis is the same contradiction. 
e Smith's result can be derived just from slope-independent visibility, so 
there may be a better way to do the ray-tracing derivation. 
e We can't construct a heightfield and just path trace it. 


e Any heightfield with a finite number of heights must violate the 
assumption of all heights being fully independent. 
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Smith masking is weird (2/2) 


e Ray-tracing derivation has odd result for downward rays (u < 0). 
e Derivations don't require u > 0, but do require u = 0. 
e Can show that A(—yp) = —A(u) - 1. 


e With 4; « à, and u > 0, we have S(éo,é1, 4) = Ga 


-A(u)-1 
. SCi, fo, =u) = (3) = ED SE, ep LL) < S(£o, T LL) 
e Visibility is asymmetric; rays traveling the same path between two heights 
are more likely to hit something going down than going up! 
e Derivation assumes rays can hit any opposing microfacet normal (V-m< 0). 
e More normals oppose downward rays than upward rays. 


e Since the probability distribution of normals is everywhere the same, the 


cumulative area of candidate normals must be greater for downward rays than for 
upward rays. 


e In other words, surface area is bigger going down than going up. 


p^ 
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Path tracing with Smith masking 


e  BRDFs must be symmetric 
To be symmetric, we need S(Z,,£s,—u) = S(£o, EL) 


e Derivation instead has S(à&, £y, —L) = PO S Fo, SW) 
1 


e Caused by A(—p) = —A(u) - 1 
To "fix", redefine A(u) so A(—p) = —A(u), without changing positive slopes: 


o e CRINES 
H Jiul 


e Conceptually, renormalize downward surface area to match upward surface area. 


For GGX, A(u) = 5 (ee) yy +0 
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Derivation that A(—u) = —A(u) — 1 


e This derivation uses the fact that P,(—q) = P;(q) 


© AC) = — f^ (q+ w)P2(q) dq 


BR 
«cs f^ Po) dq — 7 f, (a — 4 + 200 P4 (4) dq 

e 7S" aPo(q) dq — 7 f^, uPCI) dq — Mu) — 2 JP Pala) dq 
e 0- fe, Pa(q) dq — AQ) — J, Pala) dq — Juz, Pala) da 

e AW — f, Pala) dq 


e —A(u)—1 
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Handling u = 0 


e Preceding derivation assumed u = 0. If you instead assume u = 
0, you get 
P4($9) (O 
© S(T) = oot FE) Ju (4-1) P2 (2) dq) 
2 ay EPI) 
e For GGX, fj qP2(q) dq = fy ss da = 2, SO S(&, 0,7) = e "2 76v 


e This barely resembles the equation for u + 0 


_ (fo) A f (So tut) 
e S (čo. 4, T) U (2) B ( f ($0) 


e Limit as u -> 0 of the equation for u + 0 is the equation for u = 0! 
Laten) SAUE ai au) ARD 
f §o) f (šo) 
e True even though A(z) is discontinuous at 0, since 


) -7 Ju (4-H) Po(a)aq 


e Furthermore, lim, ( 
uot u07 


f §o+HT) > 1. 
0 
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Heightfield heights 


e Given D(m), it's possible to figure out the heightfield height limit. 
e We have P,(q), the 1D probability of slope q. 


2 
e GGX's cumulative probability X = ko ree = Lu 


e Solve for u to turn a random variable into a slope: u = A 


f . (2Xj—1) (2X;-1) 
e Height IS Yi pdt = YT ead = recy aera 
e All roughnesses can use same heightfield, just scaled by a. 
e Correctly says a = 0 is perfectly flat. 
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GGX heightfield probability func 


e Isummed random slopes with a = 1 to generate heightfields. 
e The random number generator is proven good with a period around 2ºº. 


2 +0.0100 | | | 
220 +0.0030 | | | | | 
224 +0.0007 TM (| | Il 
227 +0.0004 ||| | EN Il | | 
230 +0.0002 | || | | | | || | | 
me +0.0002 WERI | ha NO) VE | " 


e Height histograms were spiky with no correllstion iid runs. 
e Each run basically picked a random number of random heights to center on. 
e Uniform height distribution over +0.0002a@ seems reasonable. 


Starting to path trace GGX 


e We now have what we need to path trace GGX: 


—0.0002 < E < 0.0002 


otherwise 


e (= ju 
0 


© fŒ = fÉ P. GO dx 
é A(u) = ip w), TE SEU 


( f(Eo) y" "T 
e S(£y,u,T) = 4 M Gore 


e 2 fo) u=0 


GOC | FEB 27-MAR 3, 2017 | EXPO: MAR 1-3, 2017. #GDC17 


Intersection distance (1/2) 


e Going from height & in direction u, at what t do we hit the surface? 
e Cumulative probability of hitting the surface is 1 — S, since S is the cumulative 
probability of not hitting the surface. 
e We can pick a uniform random variable for 1 — S and solve for c 
e This is equivalent to picking a uniform random variable for S 


( f) p 3:8 
© SGo T) = 4 VED 


-r Pio) 
e 2 fo u=0 


fF (f Go) s-VAQ0)-£, 


0 
u id 


- aneaf _ 
In $75 G0 E 
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Intersection distance (2/2) 


e More convenient to use @ and d. 
e Lu -—cotÜü tT=dsin6 yur-—dcosO 
: A(8) . Va2+(1-a2) cos? 8—-|cos 6] 


2 cos 0 cos 0 #0 
(fy n 
e S(£5,0,d) = 4 Yota cos 8) 
-d a P1($0) 
e 2 f($o) cos@ = 0 
(ff Eo) S7140) — Es)tang cose #0 
did 2f ($9) o 
—In S Th GO cos6 = 0 


e For a 3D vector V, we have cos 6y = V,, making this trivial to calculate 
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Escaping rays 


e If cos8 #0, then d = (f !(f(£9) SAO) — Es)tang 
e f^! is undefined if f(z,) S714 > 1. 

e Can only happen if A(@) > 0, which is when cos6 > 0 (upward rays). 
e Fortunately, algebra shows this is when S < f(&)4®. 


e Recall that f())4 is the probability of a ray escaping when it 
starts at é& and goes in direction 8 


© If S< f(&)4™, the ray hit the viewer, not the microsurface. 
e This is our one and only path termination condition. 
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Intersection normal 


e Smith derivation assumes D(m) is independent of height. 


e So, for a vector traveling in direction T, pick any m according to 
D (m) such that T -m < 0 (ray points at surface, normal points 


away). 
e Need to pick (m, m). 
e D(m) = Es 


m(1-(1-a?) cos? Om)? 
e Since D(m) doesn't depend on azimuth m, just uniformly pick in 
[0,27). 
e Need to importance sample D(m) over hemisphere 


e k2r f? Dam) cos Om sin Om dO, = 1 
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Intersection normal 


X : 
e Y= 21 f; D(m) cos Om sin Om dom 
_ fX 2a* cos Om sin Om 
« f= Jo (1-(1-a@?) cos? 0m)? 
e Note that rom. — (1 — a?) cos? 64, = 2(1 — a?) cos Om sin 8m 


dOm 


a? I 
——3]9 (x) . . 2 
e We have the form Eres which has the solution -( : ) n 
g(x) 1-a2/ g(x) 
y = ME | x a la] 
° B 1-a?/ 1-(1-a?) cos? Om lo 7 1-a?/ \1-(1-a@2) cos? X a? 
y= ( a? ) a?-(1-(1-a?)cos? X)Ó| ^ -1*co?xX  — 1-cos? X 
° = 1-a2 a?(1-(1-a2)cos? X) J)  1-(1-a2)cos?X 1-(1-a2)cos? X 


e When X == (whole hemisphere), cos X = 0, so Y = 1 (i.e., already properly 
normalized). 
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GOC 
Intersection normal 


y = 1-cos? X 
j -~ 1- (1-02) cos? X 
Need to solve this for X given Y. 


Y —Y(1— aĉ?) cos? X = 1— cos? X 
(1—-Y(1—-2a?))co?X =1-Y 


3 1-Y 
OS^ X = ———— 
» So 1-ü-a2)Y 


= -1 1-Y 
e X — cos ( s 


m) 


E : , Z —1 
e Algebraically equivalent to Walter's result X = tan (zm 


Intersection normal 


We can now pick a random microfacet normal: 
Xo, X1 = uniform random values in [0,1] 


e dm = 27 Xo 


o A | tx 
e Om = cos — 


e m= (cosi, sin 0, , sin m sin Om , cos Om) 


e Start over if T-m 2 0. 
Have to retry because there is no closed form solution to 
importance sample GGX's D(m) directly given the constraint T - 


m « O0. 
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Intersection normal 


e Can be slow to pick a normal as T, > 1. Almost none of 
our guesses satisfy the constraint T -m < 0. 


e Dot product is cos m sin Om Ty + sin dm sin Om Ty + cos Om T, < 
0. 


e T4, COS Ọm + T, sin m < —T; cot Om 
e Set T, =rcosf and T, =r sinf 


ED + Ty = —T,*; B = some unknown angle 


e rcosfPcos¢d, +r sin fp sind, < —T, cot Om 
_ 


r tan Om mando 


the minimum value for the left is — 


e cos(B — dy) <— 
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Intersection normal 


Tz 


= == Tz 
o 1< su SO tan n > ^ 
1-Tz 
e cos?0,«1-— T, (trivial since T, acts like a sine of some angle) 
1-X š š a š 
dio m T,^ (the left is our derivation for sampling cos? 8m) 
RON 1 


© 1-Xi«(1-T/7/)-(1-7,?)-a?)x, 
© T< (1 de de a?) X, 


Ln = n 
1-(1-T?7)1-a?)  a?«T;?(1-a?) 


e X > 
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Intersection normal 


e We can now pick a random normal more efficiently: 
e Xy, X = uniform random values in [0,1] 
o dm = 27 Xo 
e If T; 5 0, shrink X, to the range [ Us 


a2 T7? (1-a?)' 


1] 


= aci | 1-X4 
e Om = cos (==) 


e m= (cos dy» sin 8n , sin $,, sin 04, , cos Om) 

e If T m« O0, return m 

e If T; » 0, negate $4. If T- m « 0 now, return this m. 
e Start over 


e This is more efficient, because at least half the values for Xy, X, generate 
valid normals given the constraint T-m < 0 
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Reflection/Transmission 


Given the microfacet normal m and incoming direction T, we can 
calculate Fresnel F to decide to reflect or transmit at the facet. 

e F=h+(1-F)(1-m-T), Schlick's famous approximation 

e We use F, = 0.02, for index of refraction = 1.33, common for dialectrics. 
Each ray starts with 1 unit of energy. If the ray's energy is above a 
threshold, we split it into reflected and transmitted parts with energies 
scaled by F and (1 — F), respectively. Otherwise, we use Russian 
Roulette to decide which path gets all the energy. 


Reflected rays continue recursively in the reflection direction. 


Transmitted rays continue in a carefully chosen random direction. 


e If the transmitted ray's energy is above a threshold, we split it into N rays 
first. 
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Transmission Direction 


e Lambertian scattering would use a cosine weighted hemisphere. 
e We tried that. The BRDF was not symmetric. 
e The problem is our effective microfacet BRDF was: 
e F «specular + (1— F) + diffuse 
e Both specular and diffuse are valid BRDFs, but F lerps between them based 
only on the incoming direction T = —L. 


e This means that swapping L and V is asymmetric in F; it replaces one vector 
with an unrelated one. 


e In short, Lambertian diffuse doesn't play nicely with a specular BRDF. 
e Fortunately, Shirley et al. solved this in 1997. 
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Transmission Direction 


e Why is Lambertian cosine weighted? 

e Lambertian scattering assumes that light enters the microsurface, bounces 
around on the inside, and then comes back out. 

e When there is a scattering event inside the surface, it assumes each 
outgoing direction is equally likely for a ray. 

e You can thus model the interior volume as having uniform beams of energy 
in every direction. The ones pointing to the surface escape. 

e But the BRDF is defined for a unit area of the surface, not the interior 
volume. A unit area on the surface cuts diagonally across the uniform beam 
exiting at an angle, so that the fraction hitting the surface is only cos@. 
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Transmission Direction 


e We had a Fresnel reflection on entering the microsurface 
volume. For symmetry, we need the same Fresnel 
reflection on exiting too. 


e Fortunately, reflection/transmission is symmetric. 


e [his means we can calculate Fresnel transmission from the 
view direction into the surface, and it is equivalent to 
calculating the Fresnel inside the surface for another vector 
that gets refracted into the view direction. 
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Transmission Direction 


e So, the probability of keeping an exiting 
direction is 1 — F(cos 6,): 

e 1— (Fo + (1— F9)(1— m-V)) 

e (1— F)(1-(1- m-V)?) 

e We need a normalization constant such that this 
integrates to 1 over all view directions. That's 
because rays reflected back into the surface will 
bounce around and get another chance to escape. 


e Since 1 — Fy is a constant, we can just absorb it as part of 
the normalization constant 
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Transmission Direction 


TU 
° k [2(1— (1 — cos @)*) sin do = 1 
e The outside the exponent is the normalization constraint for a 
BRDF 
e sin dð is the measure for integrating over the hemisphere. 


° is from integrating over the hemisphere but not depending on 
azimuth. 


e kis the normalization constant we want to find. 
e This is actually easy to solve. Just multiply out (1 — cos 6)” to 
get a polynomial in cos8. The 1’s cancel. We're left with terms 
ike: 


e a cos? 0 sing 
e Trivial integral of each term is — = cos^*! 6 
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Transmission Direction 


e The final result of the integral is: 


5 3 5 4 5 5 6 1 7 Ov 
e 2nk (—=cos 0 +5cos 0 —2cos 0 t — cos 0 — 7 cos e) 
0 


e For6,- A cos 0, = 0, and we're left with 


5 1 10 207 
2nk ($722 - 52) = ank (=) =o," = 1 
ee 
207 T 
«Interestingly, this is just 5% larger than the pure Lambertian BRDF. 
21 5 5 5 1 
e 1+ = (— cos? 0, + = 605" 0, — 2cos? 0, + ; cos? 0, — = cos" 6,) 


7 21 21 7 3 
e 1— > cos? 0, + 4 cost 0, == 008" 6, += cos 0, — o COS” 0, 


e We want to importance sample this with a [0,1] variable, so can use 1 - this. 
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Transmission Direction 


7 21 21 7 3 
e Y= > cos” 0, — = cos* 0, + = cos? 0, — ; cos" 0, + mid 0, 


e We can solve this polynomial for cos 6, to use in importance 
sampling of view directions. 


e No easy closed form solution; follow Shirley's recommendation to 
pick a good guess then improve with Newton Raphson. 


_ 0.0114813--154.4Y -13002.4Y ?  38295.9Y? 
1+1483.577+33596.4Y2+16520.2Y3 


e We then use this and 4, in [0,27] to get the outgoing 
transmission direction relative to the microfacet's normal m. 


e The first guess is y = cos0, = 
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Transmission Direction 


e The result is a combined BRDF: 


e Fx specular + (1— F) «ea —(1— m-V)*) 
e [he specular part is symmetric because it is nonzero only where m = H, which is 
where F(m -L) = F(m-V) 
e H = normalize(L + V), the half-angle vector. 
e The diffuse part is symmetric because 1-F=(1-F)(1-(1-m-L)) 
e Swapping L and V just swaps which Fresnel is entering and which is exiting, as expected. 
e Fis what fraction of rays reflect; 1— F is what fraction transmits. 
e If the ray reflects, only facets aligned to H reflect light, without loss of energy. 
e Ifthe ray transmits, it has to survive another transmission event to escape the 
surface. Once it escapes, partial absorption has tinted it by pa. 
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Transmission Observation 


e This still doesn't perfectly model subsurface effects 
e Most obviously, we ignore how Snell's law changes ray 
directions 
e Assumptions have ray directions uniformly distributed inside the surface 
e Equation assumes rays are uniformly distributed outside the surface 
e Snell's law says y, sin0, = nz sin 802, SO 0; = sin! (= sin 0.) 


e 0, approximately linear for 6, near O, then 6, changes faster as it 
approaches 7 — and 6, approaches angle of total internal reflection 


e This mes: outgoing angles are less represented near E Z than near O 


e Conveniently, that’s where transmission is weakest, so they're already less 
represented 
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Appendix Conclusion 


This now gives all the pieces needed to get the path tracing 
solution 


I picked about 64 representative zenith angles and 16 « values, 
and then shot tons of rays for each pair of settings. 


For each ray, I recursively saw what it hit and did Fresnel 
transmission/reflection based on the chosen normal. Finally, I 
bucket escaping rays into view directions. 

From this raw data, I tried tons of random equations with terms 
symmetric in L and V until I saw ones that I liked based on their 
tradeoff between computation cost and fidelity. 


